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Abstract. In an XMM-Newton raster observation of the bright local group spiral galaxy M33 we study the population of 
X-ray sources (X-ray binaries, supernova remnants, super-shells) down to a 0.5-10 keV luminosity of 10 35 erg s" 1 - more 
than a factor of 10 deeper than earlier ROSAT observations. EPIC spectra and hardness ratios are used to distinguish between 
different source classes. We confirmed the 3.45 d orbital light curve of the X-ray binary M33 X7, detected a transient 
super-soft source in M33, and searched for short term variability of the brighter sources. We characterize the diffuse X-ray 
component that is correlated with the inner disk and spiral arms. We will compare the results with other nearby galaxies. 
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1. Introduction 



The Local Group Sc galaxy M33 at a distance of 795 kpc 
(van den Bergh 1991) with relatively low inclination of 56° 
(Zaritzky, Elston & Hill 1989) and its moderate Galactic fore- 
ground absorption (Nn= 6 10 21 cm -2 , Stark et al. 1992) is 
ideally suited to study the X-ray source population and dif- 
fuse emission in a nearby spiral. The Einstein X-ray Obser- 
vatory detected diffuse emission from hot gas in M33 and 
17 unresolved sources (Long et al. 1981; Markert & Rallis 
1983; Trinchieri, Fabbiano & Peres 1988). First ROSAT HRI 
and PSPC observations revealed 57 sources and confirmed 
the detection of diffuse X-ray emission which may trace the 
spiral arms within 10' radius around the nucleus (Schulman 
& Bregman 1995; Long et al. 1996). Combining all archival 
ROSAT observations of the field, Haberl & Pietsch (2001; 
HP01 hereafter) found 184 X-ray sources within 50' radius 
around the nucleus, identified some of the sources by corre- 
lations with previous X-ray, optical and radio catalogues, and 
in addition classified sources according to their X-ray prop- 
erties. They found candidates for super-soft X-ray sources 
(SSS), X-ray binaries (XRBs), supernova remnants (SNRs), 
foreground stars and active galactic nuclei (AGN) in the back- 
ground. 



Two M33 sources are known for their outstanding X-ray 
properties (Peres et al. 1989). The brightest source (X8, lu- 
minosity of about 10 39 erg s _1 ) is the most luminous X-ray 
source in the Local Group of galaxies and coincides with 
the optical center of M33. Its time variability (Dubus et al. 
1997), its point-like nature seen by ROSAT HRI (Pietsch & 
Haberl 2000) and Chandra (Dubus & Rutledge 2002) and its 
X-ray spectrum best described by an absorbed power law plus 
diskblackbody model (e.g. Ehle, Pietsch & Haberl 2001, E01 
hereafter; La Parola et al. 2002) point at a black hole XRB. A 
possible periodicity of 106 days was not confirmed in later 
observations (see Parmar et al. 2001). The second source 
(X7) is an eclipsing XRB with a binary period of 3.45 d and 
possible 0.31 s pulsations (Schulman et al. 1993, Dubus et al. 
1997, 1999; D99 hereafter; Larson & Schulman 1997). 

Here we present XMM-Newton raster observations which 
were carried out within the telescope scientist guaranteed 
time program to survey M33 homogeneously for point 
sources with a sensitivity of 10 35 erg s _1 in the 0.5-10 keV 
band, a factor of ten deeper than previous surveys. First X- 
ray false color images combining different energy bands have 
been discussed in Pietsch (2002). 

2. Observations, data analysis and results 
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Fifteen 10 ks observations on M33 were scheduled with 
XMM-Newton (Jansen et al. 2001) between August 2000 
and January 2002. The first two observations were performed 
with the thick filter in front of the EPIC detectors (Turner et 
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Fig. 1. XMM-Newton EPIC low-background exposure map 
of the M33 raster. Contours are at and 5 ks, gray-scale steps 
every 2 ks. The maximum exposure is 35 ks in the SW M33 
disk. The optical extent of M33 is marked by the white D25 
ellipse. 



al. 2001; Striider et al. 2001) while for the rest the medium 
filter was used. The raster pointings with a spacing of ^10' 
ensure that each position within the optical D25 ellipse of 
M33 is covered at least 3 times. During the proposed 30 ks 
integration time, for each source inside the ellipse at least 
15 counts are detected in the EPIC PN camera in the 0.5-10 
keV band (taking into account the X-ray telescope and instru- 
ment response and assuming a 5 keV thermal bremsstrahlung 
spectrum). Several observations suffered from high particle 
background reducing in some cases the integration time use- 
ful for the raster imaging to below 1.3 ks. Fig. [j] shows the 
EPIC PN coverage of M33 (total integration time 1 12 ks, 92 
ks with medium filter). Several pointings have been granted 
for reobservation to homogenize the survey. 

The data analysis was performed using tools in the SAS 
v5.3.3 and FTOOLS v5. 1 software packages, the imaging ap- 
plications DS9 v2.1 and KVIEW vl.1.19, the timing analy- 
sis package XRONOS v5.19 and spectral analysis software 
XSPECvll.2. 

2.1. X-ray images 

The creation of the images of the raster survey of M33 was 
performed in several steps. We first cleaned the EPIC event 
lists of the individual observations for bad pixels and times of 
high background and calculated sky positions with respect to 
a common reference point close to the center of the galaxy. 
We then accumulated images and exposure maps for the in- 
dividual observations in the energy bands 0.2-0.5 keV (Bl), 
0.5-1.0 keV (B2), 1.0-2.0 keV (B3), 2.0-4.5 keV (B4), and 
4.5-12 keV (B5). For EPIC PN we also created "out of time 
event" (OOT) images. For EPIC PN we excluded the energy 




Fig. 2. XMM-Newton EPIC low-background image of M33 
(0.2-4.5 keV). Data of the EPIC PN, MOS 1 and MOS2 cam- 
eras have been combined (see text). White squares indicate 
ROSAT sources from the HP01 catalogue. The optical extent 
of M33 is marked by the black D25 ellipse, the extraction ar- 
eas for Fig. H as black boxes. 

range 7.2-9.2 keV from the B5 band images, as in this sub- 
band internal florescence lines lead to a spatially very in- 
homogeneous background. We smoothed all images with a 
Gaussian of 12.6" FWHM, subtracted OOT images for EPIC 
PN, added them up and corrected for exposure times. 

The 0.2^1.5 keV image of M33 including all EPIC in- 
struments (Fig. ||) shows many more sources than previously 
detected with ROSAT within the D25 ellipse and indicates un- 
resolved emission surrounding the bright source (X8) close 
to the nucleus and from the optically bright inner disk and 
southern spiral arm. 

2.2. Diffuse emission from the inner disk 

To further investigate the diffuse emission we extracted emis- 
sion profiles in rectangular boxes along the major and the mi- 
nor axes of M33 (width of 15' and 20', respectively) from 
individual energy band images of the EPIC PN (Fig. |J). The 
bright point sources can clearly be identified as peaks in the 
plots. While in the B3 and B4 band the background level be- 
yond the point sources is rather flat, there is an enhancement 
of emission in the B2 band both within ^15' from the galaxy 
center along the major axis and ~40' from the galaxy center 
along the minor axis profiles. This excess seems stronger to 
the south and to the east of the nucleus. While this in princi- 
ple may be caused by unresolved point sources that are pre- 
dominantly radiating in the soft band (like SSS, SNRs) the 
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Fig. 3. XMM-Newton EPIC PN emission profile along the 
major (width 15', top) and minor (width 20', bottom) axes 
in the (0.5-1.0; solid line), (1.0-2.0; dotted line), (2.0- 
4.5;dashed line) keV band. 



more plausible interpretation is emission from hot gas in the 
interstellar medium (ISM) of the M33 disk or halo. Such a 
component would be similar to the hot ISM measured in the 
LMC and SMC (see Sasaki, Haberl & Pietsch 2002). A de- 
tailed spatial and spectral analysis, cutting out point sources, 
is in progress. It is however clear that this will be no easy 
task due to the extended point spread function (PSF) of the 
XMM-Newton telescopes at large off-axis positions which 
heavily smear out the bright point source near the nucleus 
if observed at larger off-axis angles. The investigation of the 
diffuse emission will therefore have to concentrate on indi- 
vidual pointings where the dominant source X8 is favorably 
positioned (ie. close to the on-axis position). 



2.3. Catalogue of point sources 

In a first attempt we produced a catalogue of X-ray point 
sources in the M33 field with the source detection programs 
of the SAS. We created unsmoothed merged images, expo- 
sure maps and masks in a similar way as described in Sect. 
2.1. Background maps were calculated taking into account 
OOT events for EPIC PN. In the overlapping fields the PSF 
at a given sky position is made up by the overlay of the in- 
dividual PSFs according to the different off-axis parameters. 
The effective PSF for a source in the merged image there- 
fore certainly differs from the one that the SAS calibration 
would assume for the offset angle from the assumed refer- 
ence point of the raster. To circumvent this effect we used a 
calibration file with a PSF similar to a source at an off-axis 
angle of 6' for the total field. This is a first approximation for 
the PSF of source counts collected at a wide range of off-axis 
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Fig. 4. XMM-Newton EPIC PN 0.3-7.2 keV light curve of 
the XRB M33 X7 (integration time 500 s) folded over the 
3.45 d orbital period of D99. 

angles, and should give reasonable count rate estimates. More 
sophisticated approaches are certainly needed. 

Simultaneous detection runs included images of all en- 
ergy bands of the individual EPIC instruments. We also com- 
bined all EPIC instruments, using simultaneously 5x3 im- 
ages (five energy bands for EPIC PN, MOS1 and MOS2). 
This resulted in a catalogue of source positions, with count 
rates and hardness ratios. 

We detected more than 400 sources in the overlapping 
field of view of the EPIC instruments. While for the fainter 
sources one can determine position, flux and possibly hard- 
ness ratios (HRs), the brighter sources (> 500 counts) allow 
in addition the investigation of their spectra and time variabil- 
ity. Detailed discussion of the source catalogue and individual 
sources goes beyond the scope of this paper. In the following 
we give some first highlights. 

2.4. Spectra and time variability of bright sources 

Based on the EPIC PN data, we performed spectral fits of the 
brightest sources and produced light curves to start the source 
identification. In E01 we gave first results on the black hole 
XRB M33 X8 including an analysis of RGS spectra of this 
source. Also shown are EPIC PN spectra of the XRB M33 
X7, a SNR, and a transient short time variable SSS candidate 
that was bright in August 2000 and no longer detectable in 
July 2001 and January 2002, a behavior often seen in SSS 
(Kahabka & van den Heuvel 1997). Many sources that were 
not discussed in E01, vary in intensity between observations. 

The XRB M33 X7 was covered by eight EPIC PN fields. 
EPIC PN counts of the source region during low background 
times were integrated over 500 s and folded over the 3.45 d 
orbital period (D99). In Fig. [^variability on 500 s time scale 
as well as on the orbital period is clearly visible. The source 
intensity seems to increase towards the end of the orbit. This 
can either reflect real orbital or long term variability of X7 as 
the data in the individual phase blocks are separated by many 
orbital cycles. The time of eclipse seems to be shifted to later 
phases by 0.1, pointing at a slightly longer orbital period than 
given by D99. Better phase coverage of the X-ray eclipse will 
allow us to significantly improve on the orbital parameters. In 
a pulsation search we could not confirm the 0.31 s pulsation 
period and did not find any other periods. 
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Fig. 5. XMM-Newton EPIC PN hardness ratio plots: (top) 
HR2 against HR1, (bottom) HR3 against HR2. Shown are 
only the 75 sources with more than 100 counts and/or ROSAT 
classifications. Two candidates for XRBs are marked as filled 
dark lozenges, 2 for SSSs as filled triangles, 15 for SNRs 
as filled dark squares, 6 for foreground stars as filled bright 
squares, 3 for AGN as open squares. 



3. Source identification and classification 

In our ROSAT work on Local Group galaxies we devel- 
oped a scheme to classify detected sources by using their 
X-ray properties (mainly hardness ratios and extent, see 
HP01). With XMM-Newton we are broadening this method 
in adding information from the harder energy bands acces- 
sible with EPIC. We define hardness ratios as HRl=(b2- 
bl)/(b2+bl), HR2=(b3-b2)/(b3+b2), HR3=(b4-b3)/(b4+b3), 
HR4=(b5-b4)/(b5+b4), where "bi" is the count rate in band 
"Bi". Specifically the use of HR3 and HR4 should allow us 
to separate XRBs from AGN, in addition to the classifica- 
tion already possible with ROSAT. With the higher sensitiv- 
ity of XMM-Newton we also will be able to classify fainter 
sources. As a first attempt, Fig. |] shows X-ray "color/color" 
hardness ratio plots for the 75 sources in the EPIC PN cata- 
logue with more than 100 counts and/or with ROSAT classifi- 
cation. ROSAT SNRs, SSS, XRBs, foreground star and AGN 
candidates are in general clearly separated. However, there 
seem to be several miss-classifications within the ROSAT list 



or miss-identifications as some of the sources do not fall into 
the regions in the EPIC PN hardness ratio diagrams populated 
by these source classes. Cross-correlations with the catalogue 
of Gordon et al. (1999) yielded several new X-ray SNR and 
AGN candidates in the field. More detailed investigations are 
in progress. 

4. Conclusions 

As demonstrated above, our setup of XMM-Newton raster 
observations of M33 yields very interesting results, fully con- 
firming our expectations. We have detected more than 400 
sources down to a luminosity limit which is below 10 35 erg 
s _1 if the sources are located in M33. We showed the pos- 
sibilities to classify them from their X-ray properties alone. 
There is unresolved emission detected at energies below 1 
keV, most likely from hot gas in the inner region of the disk 
and possibly from the southern spiral arm and/or halo above, 
that needs to be followed up. 

Several of the X-ray observations were heavily affected 
by high background and will be re-scheduled to homogenize 
the survey. Identifications of X-ray sources with surveys and 
follow-up observations at other wavelengths are in progress. 
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